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ABSTRACT 

Gossiping  models  have  been  increasingly  applied  to  study  social  net¬ 
work  phenomena,  in  particular,  to  model  the  dynamics  of  social  be¬ 
havior  or  belief  through  local  interactions.  In  this  context,  this  paper 
investigates  how  the  opinions  of  social  agents  diffuse  in  a  network 
under  a  so-called  hard-interaction  model,  in  which  the  agents  inter¬ 
act  more  strongly  with  neighbors  that  share  their  beliefs  and  have 
no  influence  on  the  neighbors  whose  opinions  differ  by  more  than 
a  threshold.  We  analyze  the  convergence  properties  of  the  opinion 
dynamics  and  provide  analytical  insights  to  characterize  the  phase 
transition  from  a  society  of  radicalized  opinions  to  one  of  conver¬ 
gent  behavior. 

Index  Terms —  opinion  diffusion,  opinion  dynamics,  social  net¬ 
works,  phase  transition,  herding. 

1.  INTRODUCTION 

The  study  of  the  convergence  of  social  behavior  can  be  found  in 
many  fields.  Examples  include,  but  are  not  limited  to  the  herd¬ 
ing  behavior  [1]  in  Economics,  th q  fad  and  trend  behaviors  [2]  in 
Social  Psychology,  and  the  Bandwagon  effect  [3]  in  political  sci¬ 
ence.  A  number  of  different  approaches  to  elucidate  these  phenom¬ 
ena  have  emerged.  Two  prominent  classes  of  models  have  been  stud¬ 
ied  extensively,  namely,  Bayesian  models  and  non-Bayesian  models. 
Bayesian  models  [1,  2]  view  individuals  as  rational  agents',  opin¬ 
ions  (or  beliefs)  of  agents  are  probabilities  of  a  given  state,  condi¬ 
tioned  on  all  the  available  information;  opinions  are  updated  using 
Bayes  rule,  with  information  communicated  through  neighbors’  ac¬ 
tions.  Bayesian  models  focus  on  the  mechanism  of  propagating  in¬ 
formation  from  one  agent  to  another  for  optimal  decision-making. 
However,  because  of  the  complexity  required  to  describe  the  diffu¬ 
sion  of  opinions  in  a  network  of  rational  agents  only  simple  sequen¬ 
tial  interactions  can  be  fully  analyzed.  Models  based  on  sequential 
interactions  usually  assume  that  agents  enter  the  society/market  one 
by  one,  and  each  makes  an  irreversible  decision  by  only  observing 
the  actions  of  its  predecessors.  This  assumption  can  be  well  justi¬ 
fied  in  certain  specific  cases,  but  it  is  hard  to  generalize  to  random 
interactions  between  agents  in  a  network.  This  motivates  the  study 
of  non-Bayesian  models  [4,  5,  6,  7,  8,  9,  10]  which  use  simple  and 
heuristic  local  belief  updating  rules  to  characterize  agents’  interac¬ 
tions  in  a  network.  Non-Bayesian  models  aim  to  capture  the  opinion 
dynamics  in  a  network  and  model  how  the  initial  opinions  and  the 
underlying  social  network  structure  affect  the  alignment  of  social 
behavior.  Earlier  non-Bayesian  formulations  include  [4,  5]  which 
model  the  interactions  based  on  simple  synchronous  linear  updates: 
individuals  assign  suitable  weights  to  their  interacting  agents  on  the 
basis  of  relative  importance.  The  updating  rule  has  the  same  form  as 
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the  well  studied  average  consensus  algorithm  [11],  and  thus  the  con¬ 
vergence  analysis  is  straight  forward.  Another  class  of  non-Bayesian 
model  is  the  so  called  Hegselmann-Krause  (HK)  model  [7,  8,  9],  in 
which  agents  update  their  opinions  using  a  nonlinear  model:  inter¬ 
actions  are  performed  through  a  synchronous  update  by  averaging 
all  the  opinions  that  differ  by  less  than  a  confidence  level  e;  the  rate 
of  change  in  opinions  is  determined  by  a  convergence  parameter  p 
(that  is  constant  in  time  and  across  the  network).  Similar  studies 
have  also  explored  the  effects  of  simple  interactions  between  two 
neighboring  agents.  For  example,  Deffuant  et  al.  [10]  modeled  the 
network  on  a  square  grid,  in  which  the  agents  can  only  communi¬ 
cate  with  their  four  immediate  neighbors  and  exchange  their  beliefs 
with  a  fixed  weight  parameter  if  the  distance  in  opinion  is  smaller 
than  a  given  threshold.  Weisbuch  [12]  extended  this  simple  lattice 
topology  to  a  scale  free  network  model  and  also  examined  the  con¬ 
vergence  property  using  a  heterogeneous  constant  threshold.  The 
focus  of  these  papers  is  on  the  modeling  of  social  interactions  and 
on  how  to  reach  opinion  consensus.  Analysis  of  the  above  nonlinear 
models  is  generally  carried  out  by  extensive  computer  simulations, 
but  explicit  mathematical  results  are  limited. 

In  this  paper,  we  explore  a  theoretical  framework  called  the 
hard-interaction  model,  to  provide  analytical  insights  on  the  asymp¬ 
totic  behavior  of  a  social  group  in  relation  to  the  agents’  initial 
opinion  profile  and  the  underlying  network  structure.  The  hard- 
interaction  model  generalizes  the  binary  decision  making  (decision 
between  two  alternatives)  as  assumed  in  [7,  8,  9,  12]  to  the  multi¬ 
alternative  decision  making  (decision  between  multiple  alternatives). 
Specifically,  rather  than  restricting  the  opinions  to  lie  in  a  bounded 
(real)  interval,  we  treat  each  agent’s  opinion  as  a  vector  of  probabil¬ 
ities;  each  element  of  the  opinion  vector  represents  the  probability 
that  a  certain  alternative  is  true.  Furthermore,  we  extend  the  HK 
model  by  introducing  a  trust  function  p  similar  in  spirit  to  the  con¬ 
vergence  parameter  p  defined  in  [7,  8,  9,  10,  12],  but  allowing  the 
trust  function  p  to  vary  with  the  opinion  distance  between  the  in¬ 
teracting  agents.  Hence,  p  is  time  varying  and  its  value  depends  on 
how  the  distance  is  defined.  Finally,  we  let  the  underlying  social 
network  to  be  any  arbitrarily  connected  network.  Agent  interactions 
are  pairwise  random  encounters.  Opinion  updates  are  as  follow:  the 
distance  between  the  opinions  of  two  interacting  agents  decreases 
only  if  the  previous  distance  between  their  opinions  is  smaller  than 
a  threshold  r,  otherwise  the  opinions  remain  unchanged;  the  rate 
of  change  of  the  opinion  distance  is  governed  by  the  time-varying 
function  p.  (Note  that  if  p  is  kept  constant  and  Li-norm  is  used  to 
measure  the  opinion  distance,  then  our  model  is  analogous  to  the  HK 
model  for  binary  decision  making.)  The  contribution  of  this  paper 
is  three-fold:  (i)  we  provide  analytical  insights  on  the  asymptotic 
behavior  of  the  social  group  and  show  the  existence  of  a  phase  tran¬ 
sition  from  diverging  clusters  of  opinions  to  herding;  (ii)  we  prove 
that  a  necessary  condition  for  society  to  herd  with  probability  1  is 
that  the  threshold  r  be  strictly  greater  than  the  expected  initial  opin- 
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ion  distance,  averaged  over  the  edge  probability  distribution;  (iii) 
we  also  show  that  if  the  social  fabric,  represented  by  the  probability 
distribution  of  pairwise  communications,  is  random  and  independent 
of  the  initial  belief  distribution,  then  the  social  fabric  does  not  effect 
the  phase  transition  towards  herding.  All  our  findings  are  validated 
numerically. 

2.  INTERACTION  MODEL 

In  our  model  there  are  V  =  {1,  2,  •  •  •  ,  n}  social  agents.  They  can 
interact  with  other  agents  at  random,  but  only  if  they  are  connected 
through  an  edge  of  a  fixed  communication  graph  Qc  —  (V,  £c)» 
where  8C  is  the  set  of  edges.  Throughout  the  paper,  Qc  is  assumed  to 
be  connected  and  Nj  denotes  the  set  of  neighbors  of  agent  j,  whose 
cardinality  \J\fj\  is  called  the  degree  of  agent  j.  We  model  the  j th 
agent’s  opinion  as  a  ^-dimensional  vector  Xj  =  [xj i,  •  •  •  ,  Xjq\  in 
X  =  {x  =  [xi,--  -  ,xq]  |  xi  =  landau  G  [0,1]}.  Al¬ 

though  this  is  not  critical  for  the  technical  derivations,  we  can  in¬ 
terpret  Xk  as  a  conditional  probability  P(Qk\€i),  where  6 1,  •  •  •  ,0q 
are  q  <  n  possible  outcomes  of  an  experiment  0,  and  £$  is  the  ith 
agent’s  private  information  on  0,  generated  from  a  certain  probabil¬ 
ity  measure  P(£|0).  The  vector  Xj  [0]  is  the  belief  of  agent  j  prior  to 
any  interaction,  while  Xj  [k]  represents  the  belief  after  k  interactions. 
In  our  simulations,  the  initial  beliefs  Xj  [0]  of  the  agents  are  drawn 
from  a  uniform  distribution  over  X. 

We  introduce  a  proper  distance  function  d(xi,  Xj)  :  X  x  X  — * 
M+  (R+  is  the  set  of  real  non-negative  numbers)  to  measure  the 
degree  of  agreement  between  agents  i  and  j.  With  respect  to  the 
norm  \\xj\\  :=  d(xj,Xj ),  we  assume  that  the  set  X  is  bounded, 
i. e. ,  sup  ^  \\xj\\  <  oo\/xj  G  X  and  Vj  G  V.  Hence,  d(xi,Xj)  < 
2supi||aJi||  :=  dmax  using  the  triangular  inequality.  Agents  in  the 
network  interact  at  random  and  change  their  beliefs.  We  define  a 
time-invariant  vector  p  whose  ith  element  pi  is  the  probability  of 
node  i  initiating  an  interaction  and  the  stochastic  matrix  P  whose 
(i,j) th  element  Pij  denotes  the  probability  that  node  i  will  choose 
to  interact  with  node  j.  Matrix  P  is  assumed  to  have  the  same  struc¬ 
ture  as  Qc :  if  an  edge  (i,  j )  does  not  exist  in  Qc,  then  Pij  =  0.  We 
assume  Pa  =  0  and  that  the  matrix  P  is  fixed.  Let  x  and  d'  denote 
the  variables  after  a  generic  update  and  dij  [k]  denotes  the  distance 
d(xi[k\,Xj[k])  after  k  network- wise  interactions  have  occurred.  It 
is  assumed  that  the  distance  between  beliefs  after  an  update  cannot 
be  larger  than  that  before  the  update.  Formally,  the  degree  of  change 
in  belief  is  captured  by  a  discontinuous  function  p(d)  according  to 
the  following  nonlinear  model: 

(al)  dij[k  +  1]  =  (1  ~ekp(dij[h]))dij[k]  ,  (1) 

where  €k  denotes  the  step-size.  For  technical  reasons,  we  need 
(a2)  ek  :  ek  ->•  oo  ,  ek  <  oo  • 

One  example  of  an  interaction  that  leads  to  (1)  is  d(x'i,x'j)  = 
d(xi ,  Xj)  —  d(xi ,  x'i)  —  d(xj ,  Xj),  which  geometrically  means  that 
the  beliefs  are  moving  closer  through  the  shortest  path  connecting 
them.  Regarding  (1),  the  following  assumptions  hold. 

(a3)  p(d)  is  a  non-increasing  function  of  d. 

(a4)  p{d)  is  C2 -differentiable  and  0  <  p{d)  <  l/ek\/d  G  [0,  r). 
(a5)  r  :  Vd  >  r  — >  p(d)  =  0  . 

(a6)  p(0)/p(r~)  <  ft  <oo  . 

(aT)  h(d)  =  p(d)d  is  concave  Vd  G  [0,  r). 

Specifically,  p(d)  represents  the  degree  of  change  in  opinion  dis¬ 
tance  after  each  interaction.  (a3)  implies  that  agents  interact 
strongly  (/.  e. ,  large  displacement  of  belief)  with  neighbors  that 
share  their  beliefs,  and  (a5)  indicates  that  agents  have  no  influence 
on  the  neighbors  whose  opinions  differ  by  more  than  r,  where  r  is 


a  measure  of  how  openminded  a  society  is.  For  the  communication 
model,  we  have  Pij  >  0  (fj)  G  Sc.  Define  P  to  be  an  n  x  n 
matrix  containing  the  probability  that  the  pair  (i,j)  performs  an 
exchange,  i.e.,  Pij  =  PiPij  +  PjPji  •  The  uniform  communication 
model  corresponds  to  the  selection  of  homogenous  rates  pi  —  1  /n, 
and  Pij  =  1  /  \Mi  | ,  uniform  across  neighbors  Wk. 

We  say  that  the  network  attains  consensus  ( herding )  iff  Vz,  j  G 
V,  d(xi[k\,Xj[k\)  =  0  for  some  k.  This  does  not  imply,  how¬ 
ever,  that  all  agents  are  certain  about  a  specific  outcome  of  an  ex¬ 
periment  0,  i.e.,  lim/c^oo  Xi^[k\  —  xf  S(£  —  £*)  for  some 
l*  G  [1,  •  •  •  ,  q\.  In  other  words,  consensus  is  achieved  when  all 
agents  have  the  same  belief  vector  w.r.t  the  distance  measure,  but 
need  not  believe  in  only  one  outcome.  If  network-wide  consen¬ 
sus  is  not  reached,  it  is  possible  that  the  network  evolves  into  non¬ 
interacting  sub-networks  (herds),  each  of  which  is  internally  in  con¬ 
sensus;  we  call  such  a  process  radicalization. 

3.  ANALYSIS 

Under  (a2),  using  stochastic  approximation  theory  [13],  we  can  map 
(1)  onto  the  ordinary  differential  equation  (ODE), 

dij  =  -p(dij)dij  ,  (2) 

in  which  dij  is  the  derivative  of  dij  with  respect  to  a  continuous 
variable  t  replacing  the  discrete  index  k.  For  the  sake  of  notational 
convenience,  we  do  not  explicitly  show  that  time  variable  t  in  (2)  and 
the  rest  of  this  section.  Let  d  be  the  average  of  dij  over  the  edges 
(i,j)  €  Sc,  i.e.  d  :=  J2(i,j)eEcpijdij •  Then’  usin£  (2)>  we  £et 

d  =  -V  Pijp(dij)dij  .  (3) 

yT'iJ  Jci-c 

Let  S  —  {(i,j)  G  Sc\dij  <  d}.  Under  (a3)  and  V(z,j)  G  S, 
p(dij)  >  p(d)  and  thus  d  <  —  p(5)  which 

equals  to  zero  when  either  p(d)  =  0  or  jj^sPijdij  =  0.  If 
both  are  positive,  then  the  system  has  to  converge.  However,  as  S 
changes  dynamically,  deriving  a  sufficient  condition  is  not  trivial:  in 
remark  1,  we  show  through  a  counterexample  that  p(d)  >  0  is  not 
a  sufficient  condition  for  convergence.  We  next  establish  a  lower 

bound  on  the  rate  change  d. 

Property  1  Under  (a 3),  (ah)  —  (a7),  when  d  <  r,  the  system  in  (3) 
is  lower  bounded  by  d>  —f3p(d)d  . 

Proof  See  Appendix.  ■ 

Based  on  the  above  result,  it  is  expected  that  (3)  will  not  con¬ 
verge  if  the  lower  bound  does  not  converge.  For  convenience,  we 
use  b  instead  of  d  to  denote  the  distance  for  the  lower  bound  system 
whose  dynamic  is  expressed  as  b  =  —/3p(b)b.  In  fact,  the  dynamic 
of  b  locally  resembles  the  form  of  the  logistic  equation  which  is  in¬ 
vestigated  in  the  next  lemma. 

Lemma  1  Assume  that  the  network  is  connected.  Under  (al)  — 
(a7),  the  system  b  =  —/3p(b)b  converges  ifr>  &(0). 

Proof  See  Appendix.  ■ 

Given  that  (3)  is  lower  bounded  by  —f3p(d)d  and  Gc  is  con¬ 
nected,  we  use  the  previous  lemma  to  establish  a  necessary  condition 
for  (3)  to  converge  under  the  proposed  model  in  (1). 

Lemma  2  Assume  that  the  network  is  connected.  Under  (al)  — 
(a7),  a  necessary  condition  for  the  system  in  (3)  to  converge  almost 
surely  is  r  >  d[0]. 
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Proof  From  Lemma  1,  the  system  b  =  —/3p(b)b  converges  if  r  > 
6(0).  Because  (3)  is  lower  bounded  by  this  system,  given  a  threshold 
r,  (3)  will  not  converge  if  r  <  d(0),  or  d[0]  w.r.t.  the  discrete-time 
system  in  (1).  ■ 

Lemma  2  indicates  that  the  system  will  converge  if  the  threshold 
is  above  d[0],  the  average  initial  distance  between  agent  pairs  that 
can  interact.  In  other  words,  if  the  system  is  open-minded  enough 
(relative  to  this  initial  dissonance  in  opinions),  then  the  system  will 
converge.  Interestingly,  what  we  observe  from  the  numerical  results 
(shown  in  the  next  section)  is  that  the  interaction  model  exhibits  a 
phase  transition  from  radicalization  to  herding  whenever  r  is  chosen 
sufficiently  above  d[0].  The  following  remark  is  in  order: 

Remark  1  Let  there  be  two  groups  H±  and  H2,  where  agents  have 
opinions  with  zero  distance  within  each  group,  but  V(i,  j)  :  i  E 
H^j  €  III  ,  dij  =_t  +  e.  Thenif'Eii,j)eHlxH2Pi3  < 
it  is  easy  to  see  that  d[ 0]  <  r.  This  network  cannot  converge  be¬ 
cause  the  two  groups  do  not  communicate,  as  the  rate  p(dij )  = 
0  V(i,j)  :  i  G  H\,j  E  H2.  However,  p  (d[0])  >  0  which  proves 
that  it  is  not  a  sufficient  condition  for  herding.  Note  that  in  this  case, 
=  0- 

4.  HOW  COMMUNICATION  RATES  AFFECT  HERDING 

In  the  following  lemma,  we  show  that  social  fabrics  will,  on  average, 
exhibit  the  same  phase  transition. 

Lemma  3  Let  dc  [0]  be  the  expected  initial  opinion  distance  between 
any  two  agents  in  the  set  V.  If  the  connected  graph  Gc  and  Pij 
are  random,  and  they  are  independent  of  the  initial  opinion  distance 
dij[ 0],  then  E  [d[ 0]]  =  dc[0]  ,  with  respect  to  Gc  and  the  initial 
distance  distribution. 

Proof  See  Appendix.  ■ 


and  Qc  is  fully  connected.  (Note  that  Lemma  1  and  2  hold  for  any 
connected  network.  The  choice  of  using  a  fully  connected  network 
is  arbitrary.)  The  plot  is  generated  by  averaging  over  300  runs.  Each 
starts  with  an  uniformly  distributed  random  initial  beliefs  and  ends 
after  2000  random  interactions  using  the  uniform  communication 
rate  {i.e.,  pi  =  1/n  and  Pij  —  l/(n  —  1))  and  a  function  p(d)  which 
equals  to  1  for  0  <  d  <  r  and  zero  otherwise.  Note  that  f3  —  1 
in  this  case.  From  Fig.  1,  one  observes  that  there  exists  a  phase 
transition  from  a  society  of  radicalized  opinions  to  a  society  with 
consistent  opinion  at  a  critical  threshold.  The  network  converges 
with  probability  one  if  r  is  above  a  threshold  which  is  observed  to 
be  around  0.75.  We  have  learned,  from  Lemma  2,  that  a  necessary 
condition  for  the  system  to  converge  is  that  r  >  d[ 0],  where  d[0] 
is  the  expected  opinion  distance  at  k  =  0  and  is  represented  by  the 
dotted  line.  For  all  r  <  d[0],  with  a  high  probability,  the  system  will 
not  converge  as  n  increases. 


Fig.  1.  Phase  transition  (top)  of  a  society  from  radicalized  beliefs  to 
a  consistent  belief.  Histogram  (bottom)  of  300  asymptotic  opinion 
profiles  with  n  =  100  at  r  —  0.4,  0.6,  0.8. 


In  our  simulations,  we  will  verify  this  analytical  result  and  compare 
it  with  a  heuristic  that  seeks  to  decrease  the  initial  opinion  distance 
d[0].  Specifically,  we  re-distribute  the  communication  rates  by  uni¬ 
formly  assigning  positive  interaction  rate  to  the  neighboring  agents 
(i,j)  G  Sc  who  have  similar  initial  beliefs  (below  the  threshold)  and 
zero  interaction  rate  to  the  agents  whose  initial  beliefs  are  above  the 
threshold.  Formally,  denote  by  Mi  the  neighbors  of  agent  i  in  the 
set  £eff[0]  =  G  Sc | dij  [0]  <  r}  and  \M%  \  its  cardinality.  The 

uniform  communication  over  £eff  [0]  is  defined  as 


f  t4tt  if(*>i)  e£eft[0] 

\  0  otherwise 


(4) 


This  redistribution  causes  the  interaction  probability  to  be  correlated 
with  the  opinion  distance  resulting  in  a  twofold  benefit:  it  decreases 
the  expected  initial  distance  d[0]  which  is  closely  related  to  the  crit¬ 
ical  threshold,  and  the  communications  are  not  wasted  on  the  agents 
that  have  little  influence  on  each  other,  if  at  all. 


5.  SIMULATIONS 


5.2.  Rate  of  Interaction 

Fig.  2  compares  the  phase  transitions  for  the  original  uniform  rate  of 
interaction  in  Sc  and  the  heuristic  described  in  (4)  over  three  random 
geometric  networks  (RGG)  of  n  =  50  agents.  (Note  that  the  choice 
of  RGG  is  arbitrary.)  Each  has  a  distinct  radius  of  communication 
a  G  [0,1].  In  particular,  a  =  1  represents  a  fully  connected  Qc 
and  a  =  0  corresponds  a  fully  disconnected  network.  Distance  is 
measured  in  the  L2-norm.  The  plot  is  generated  by  averaging  over 
600  runs  with  uniformly  distributed  random  initial  beliefs  over  A7  of 
dimension  3  after  k  =  2000  random  interactions.  Observe  that  the 
performances  for  a  =  1,0.8  and  0.6  are  similar  when  the  underlying 
network  is  connected  and  independent  of  dij  [0] .  It  is  consistent  with 
our  findings,  given  in  Lemma  2  and  Lemma  3.  In  contrast,  when 
Pij  is  correlated  with  dij  [0]  through  the  rate  re-distribution  (4),  the 
probability  of  forming  a  convergent  belief  increases.  For  example, 
with  r  =  0.5  in  a  fully  connected  network,  i.e.,  a  =  1,  the  probabil¬ 
ity  of  the  agents  forming  a  convergent  belief  increases  from  0.15  to 
0.95. 


5.1.  Phase  Transition 

Fig.  1  shows  the  asymptotic  algebraic  connectivity  of  a  graph 
Qeft[k]  :=  {V,£eff[fc]}  where  Ses[k]  =  {( i,j )  G  Sc\dij[k\  <  r}, 
w.r.t.  r  for  n  =  50,  n  =  100  and  n  =  200  after  k  =  2000 
network- wise  interactions.  Distance  is  measured  in  the  L2-norm 


6.  CONCLUSION 

We  proposed  a  hard-interaction  model  for  social  networks  and  ana¬ 
lyzed  its  convergence  properties  in  terms  of  the  expected  initial  opin¬ 
ion  distance.  By  describing  the  opinion  dynamics  as  the  change  in 
distance  between  opinions,  we  demonstrated  the  existence  of  a  phase 
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Fig.  2.  Performance  comparison  between  uniform  random  and  the 
heuristic  communication-rate  schemes. 


transition  w.r.t.  the  opinion  threshold  and  proved  that  the  critical 
threshold  is  lower  bounded  by  the  expected  initial  opinion  distance. 


7.  APPENDIX 


1  —  W)))  >  i  anci  (5)  converges.  The  exponential  rate  of 

convergence  equals  R(t).  Moreover,  since  6(f)  is  a  monotonically 
decreasing  function  of  f ,  an  equivalent  condition  for  convergence  is 
6(0)  <  r.  On  the  contrary,  when  6(f )  =  r,  we  get  K(t)  =  1  be¬ 
cause  p(r)  =  0.  For  6(f)  >  r,  both  p  ( 6(f ))  and  p  ( 6(f ))  are  zero, 
which  implies  R(t)  =  0  in  (5)  and  bifurcation  occurs.  Thus  the 
system  may  not  converge  when  b(t)  >  r  for  Vf . 

Lemma  3:  Define  dc  [0]  :=  E {dij  [0] }  to  be  the  expected  initial 
distance  between  any  agents  in  V.  Since  the  rates  Pij  are  random 
and  independent  of  dij  [0] ,  then 

E  [d(  0)]  =  E  [E  {iJ)eEcPi^  [dij  [0]  I  Pij] } 

=  E  |^E(Ii;)cecE’,.j]  d,  [0]  =  dc [0]  . 


Property  1:  Let  £eff  =  {(hj)  G  Sc\dij  <  r}  and  £°ff  = 
G  Sc\dij  >  r}  be  the  complement  of  £eff  in  £c-  De¬ 
fine  Peff  :=  an<^  ^eff  5^(i,j)e£eff  Pijdij' 

Under  (a7),  Jensen’s  inequality  and  the  relation  des  <  d  yield 

P^h(dii)  -  Smce 

p(dim)  =  0  for  V(^ra)  G  £eff>  hence  = 

J2(ij)e£effpijHdij)  <  dp(de ff)  .  Since  deff  <  d  <  r,  un¬ 
der  (a3)  and  (a6),  p(de&)[p(d)  <  p(0)/p(r_)  <  /3.  Hence, 
T.(i,j)££cPijKdij)  <  Pp(d)d. 


Fig.  3.  Vector  Field:  the  point  $(s;  t)  moves  toward  the  stable  fixed 
point  £i  =  0  if  ^(0;  t)  <  K.  Otherwise,  it  moves  away  from  £i. 

Lemma  1:  Suppose  that  b(t  +  s)  is  in  a  neighborhood  of 
b(t)  provided  that  s  is  small.  Hence,  Taylor’s  expansion  gives 
p  (b(t  +  s))  «  p  ( b(t ))  +  {b{t  +  s)  —  b(t))  p  (6(f)).  Then  system 
b(t  +  s)  =  —/3p(b(t  +  s))b(t  +  s )  becomes 


b(t  T  s 


-fip(b(t))b(t  +  s ) 


if  p(b(t))  =  0; 
if  p  (6(f))  <  0. 


where  f(s;  t)  =  /3  [p  (b(t))  -  b(t)p  ( b(t ))  (l  -  bit  +  s) 

and  p(6)  :=  dp/db.  In  the  first  case,  p(6(f))  is  locally  constant 
and  hence,  the  local  rate  of  convergence  around  6(f)  is  exponential 
and  is  equal  to  p  ( 6(f ))  when  p  (6(f))  >  0,  i.e.,  b(t)  <  r.  For  the 
second  case  when  p(6(f))  <  0,  define  ^(s;f)  :=  6(f  +  s)/6(f). 
With  respect  to  s,  the  dynamics  of  £  become 


£0;f)  = 

"b(t)b{t  +  s)=  PRms;t)(i 

,  (5) 

where 

R(t)  =  p(b(t)) -b(t)p{b{t))  , 

(6) 

X(t)  —  1  +  ^  (kW)  >  2 

[)  1+ ~b(t)p(b(t))>1  ■ 

(7) 

Note  that  the  dynamics  of  £(s;f)  resemble  the  logistic  equation. 
There  exists  two  equilibria  at  £i  =0  (stable)  and  £2  =  K(t)  (un¬ 
stable).  When  R(t)  >  0,  as  shown  in  Fig.  3,  the  system  will  con¬ 
verge  if  K(t)  >  £(0;f)  =  1.  Indeed,  if  6(f)  <  r,  then  K(t)  = 
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